CERN-TH.6625/92 
BUTP-92/38 



RADIATIVE SEMILEPTONIC KAON 

DECAYS 

J. Bijnens^ and G. Ecker" 

CERN 
CH-1211 Geneva 23 

and 

J. Gasser* 

Institut fiir Theoretische Physik 
Universitat Bern, Sidlerstrasse 5 
CH-3012 Bern 



Abstract 

We evaluate the matrix elements for the radiative kaon decays — > l^i^ij, 
and K —>■ nlui'-f {I, I' = e,fi) to next-to-leading order in chiral pertur- 
bation theory. We calculate total rates and rates with several kinematical cuts 
and confront the results with existing data. Measurements at future kaon facilities 
will allow for a more detailed comparison between theory and experiment. 



t Address after 1 November 1992: NORDITA, Blegdamsvej 17, DK-2100 Copen- 
hagen, 
Denmark. 

^ Permanent address: Institut fiir Theoretische Physik, Universitat Wien, Boltz- 

manngasse 5, A-1090 Wien, Austria 
* Research supported in part by Schweizerischer Nationalfonds 



CERN-TH.6625/92 
August 1992 



1 Introduction 



There are two main reasons for undertaking a systematic study of radiative 
semileptonic kaon decays at this time: 

• Very intense kaon beams will soon become available, which will increase 
the existing, rather meager, data sample dramatically. In fact, the present 
investigation was triggered by the experimental program foreseen at the $ 
factory DAFNE at present under construction in Frascati 

• The standard model can be tested in those processes to next-to-leading 
order in the chiral expansion without any further assumptions. All low- 
energy constants involved have already been determined phenomenologi- 
cally, allowing for unambiguous predictions for a large number of available 
channels. As is the case in general, there is of course no guarantee that 
effects of OijP) or higher in chiral perturbation theory (CHPT) 0, |], Q 
are small. However, the full amplitudes of 0{p^) for radiative semileptonic 
kaon decays can be directly confronted with experiment. 

The amplitudes at leading order in CHPT coincide with the current alge- 
bra amplitudes of pre-QCD times. At next-to-leading order, the amphtudes of 
O(p^) carry additional dynamical information. The one- loop amplitudes entering 
at 0{p^) are in general divergent and must be supplemented by local countert- 
erm amplitudes depending on a number of (renormalized) low-energy coupling 
constants @]. The last part is due to the chiral anomaly which contributes 
to all the i^'-decay channels considered here. At the mesonic level, all anomalous 
amplitudes can be derived from the Wess-Zumino-Witten functional P] with- 
out any free parameters. Radiative semileptonic decays thus offer a number of 
possibilities to experimentally investigate one of the fundamental aspects of the 
standard model. 

The present study completes earlier work by Donoghue and Holstein [0, § 
who did not calculate the loop amplitudes. Although non-leading in I/Nq, where 
Nc is the number of colours, the importance of loop contributions depends very 
much on the specific process under consideration. In fact, radiative semileptonic 
transitions offer a wide range of possibilities in this respect: 

• Except for the renormalization of the kaon-decay constant, there is no loop 
amplitude for Ki2'y. 

• For Ki2i'+i'', the one-loop contribution is contained in the electromagnetic 
form factor of the charged kaon. This form factor is known to be dominated 
by the coupling constant Lg at low energies 0. 
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The loop amplitude is sizeable for Ki^,^ and interferes destructively with 
the local counterterm amplitude. It is essential for a meaningful comparison 
with experiment. 



The outline of the paper is as follows. In Sect. || we recapitulate CHPT to 
next-to-leading order to the extent needed in subsequent sections. The decay 
amplitudes for Ki2^ are calculated in Sect. 3 to 0{p'^) in CHPT. The decays 
Ki2i'+i'-, T^i2e+e- and Ki^^ are treated in Sects. 4 and 5, respectively. Conclusions 
are drawn in Sect. 6. Appendix A contains a summary of our notation and of 
integrals appearing in the various one-loop amplitudes. The decomposition of the 
tensor amplitude for Ki2r^ into invariant amplitudes is discussed in Appendix B. 
Finally, Appendix |C| contains a derivation of Ward identities without using formal 
manipulations of T-products and current commutators. 



2 CHPT to next-to-leading order 

Chiral perturbation theory is a systematic approach to formulate the standard 
model as a quantum field theory at the hadronic level. It is characterized by an 
effective chiral Lagrangian in terms of pseudoscalar meson fields (and possibly 
other low-lying hadronic states) giving rise to a systematic low-energy expansion 
of amplitudes g |, 

In the formulation of Ref. [Q, one considers the generating functional 
Z[t>,a, of connected Green functions of quark currents associated with the 
fundamental Lagrangian 

^ = ^QCD + + -f^a^)q - q{s - i-f5p)q . (2.1) 

Here, CPqcd is the QCD Lagrangian with the masses of the three light quarks 
set to zero. The external fields f^, a^, s and p are Hermitian 3x3 matrices 
in flavour space. To describe electromagnetic and semileptonic interactions, the 
relevant external gauge fields of the standard model are Q 



r 



v^ + a^ = -eQA^ (2.2) 



l, = v,-a, = -eQA^--=-^— (H^;r+ + h.c.) 



Q = ^diag(2, -1, -1), T+ = I 



Vud Vus 






"'^We adopt the present conventions of the Particle Data Group |10 
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where the Vij are Kobayashi-Maskawa matrix elements. The quark-mass matrix 



M = diag(m„, m^, 



(2.3) 



is contained in the scalar field s(x) 



The generating functional Z admits an expansion in powers of external mo- 
menta and quark masses (CHPT). In the meson sector at leading order in CHPT, 
it is given by the classical action 



^2 = 1 d'^xC2iU,v,a,s,p) , 



(2.4) 



with £2 the non-linear a model Lagrangian coupled to the external fields v, a, s, p: 



where 



D^U = d^U -ir^U + iUl^, X = 2Bo{s + ip), 
and (A) stands for the trace of the matrix A; U is a unitary 3x3 matrix 

f/^f/ = l, det[/ = l, 



(2.5) 
(2.6) 



which incorporates the fields of the eight pseudoscalar Goldstone bosons. A con- 
venient parametrization of it isQ 



V2 Vg 



U = exp 



-TT 



-K+ 



\ 



V 



71 

K- 



V2 

V6 ; 



(2.7) 



The functional Z2 is invariant under local SU{?))l x SU{?))r transformations: 
Z2{v' ,a\s' ,p') = Z2{v,a,s,p) , (2.8) 

with 



gdv^ - a^)gt + igLd^,gt 
gR{s + ip)gt ■ 



s + ip 

The matrix U transforms as f/ — > gnUg^ under ( p.9| 



(2.9) 



^We follow the Condon-Shortley-de Swart phase conventions. 
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The parameters F and Bq are the only free constants at 0{p^): F is the pion 
decay constant in the chiral hmit [cf. Eq. (|2.17|) ] : 



F^ = F[l + 0{m,^ark)] = 93.2 MeV , (2.10) 

whereas Bq is related to the quark condensate: 

{0\uu\0) = -F^Bo[l + Oim,^ark)] ; (2.11) 

Bq always appears multiplied by quark masses. At O(p^), the product Boniq can 
be expressed in terms of meson masses because 

= 5o(m„ + ma) + 0{M^) . (2.12) 

The Lagrangian (|2.5| ) is referred to as the effective chiral Lagrangian of 0{p^). 
The chiral counting rules are the following: the field U is of 0(p°), the derivative 
dfj, and the external gauge fields v^, are terms of 0{p), and the fields s,p count 
as 0{p^). 

At order the generating functional consists of three different classes of 
contributions |Q : 

i) The one-loop diagrams generated by the lowest-order Lagrangian (|2 

ii) An explicit local action of order p^. 

iii) A contribution to account for the chiral anomaly. 

The contributions from categories (i) and (ii) are invariant under 
whereas those from (iii) are not. Below we will use the Bardeen form |^ of this 
non- invariant part. In this convention, the generating functional is still invariant 
under local gauge transformations generated by the vector currents, 

Z{v\a' ,s\p') = Z{v,a,s,p) ;gR = gL- (2.13) 
The local action of 0{p^) [class (ii)] is generated by the Lagrangian £4 0]: 

£4 = ...+ U{D^U^D'^U){x^U + xU^)+ L,{D^U^D^^Uix^U + U^x)) 

+ ■■■- tL.iF^'D.UDM^ + F^D^U^DM) + L,^{U^F^^'UFl,u) , 

(2.14) 

where 

F^" = d'^r'' -d^r" -i[r^',r''] (2.15) 
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Table 1: Phenomenological values for the coupling constants Ll{Mp) [i = 
4,5,9,10). The quantities Fj determine the scale dependence of the L^fi) ac- 
cording to Eq. ( |2.16| ). 



i 


mMp) X 10^ 


Ti 


4 


-0.3 ± 0.5 


1/8 


5 


1.4 ± 0.5 


3/8 


9 


6.9 ± 0.7 


1/4 


10 


-5.5 ± 0.7 


-1/4 



and where we have only written down the terms necessary for our purposes. The 
low-energy couphngs L4^, L^, Lg, Liq are divergent. They absorb the divergences 
of the one-loop graphs leading to amplitudes depending on renormalized, finite 
couplings LKfi) with the following scale dependence 

^{^^2)=^{^^,) + -^^\n^ . (2.16) 
167r^ fi2 

Observable quantities are independent of the scale fi, once the loop contributions 
are included. The coefficients Fj are displayed in Table |l| together with the phe- 
nomenological values of the renormalized coupling constants L[(Mp). In fact, the 
constants L4 and L5 will never appear explicitly in our amplitudes because we 
will use the relations 

Ml ^ Mp Fk , , 

which are valid at order in the isospin symmetry limit m„ = m^. 

We now turn to point (iii) above. The contributions of the chiral anomaly 
to CHPT amplitudes of 0(j9^) can be derived from the Wess-Zumino-Witten 
functional . Here, we will only write down the pieces relevant for the radiative 
decays under consideration. These terms can be expressed as two anomalous 
Lagrangians of 0{p^) 

CanomiW^) = - -^J^-—e^'"'^W;^ F^Jt^ {u^ D ^U, Q + ^U^ Qu]) 
4v2vr smOw \ I Z J / 

+ h.c. (2.18) 
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Canomi^h) = - j^E^"''' A, {Q [d.UdM^ d,UU^ - d,U^ d^dpU^ u]) 



e'^'P'^A^iQd^^d.^dp^) + ... , (2.19) 



where -D^t/ = df^U + ieA^[Q, U] is the covariant derivative with respect to elec- 
tromagnetism only and F^^ = d^A^, — d^A^^^ is the electromagnetic field strength 
tensor. 

We are now in a position to calculate the amplitudes for radiative semilep- 
tonic K and vr decays to O(p^) in CHPT. We will actually calculate the various 
amplitudes for an external W, which then turns into a lepton and its associ- 
ated neutrino. For the low energies that are relevant here, we can neglect the 
momentum dependence of the W propagator. 



3 Radiative K12 decays 

We consider the Ki2-y decay 

K+ip)-.i+{piMpMq) [K12,] , (3.1) 

where / stands for e or /x, and 7 is a real photon with = 0. Processes where 
the (virtual) photon converts into a e+e~ or pair are considered in the next 

section. The K~ mode is obtained from (|3.1| ) by charge conjugation. 



3.1 Matrix elements and kinematics 

The matrix element for — > l^ui'y has the structure 

T = -iGFeV:,e; {FkL^ - H^^'^Q (3.2) 

with 

= {p-qr = {pi+Pur. (3.3) 

Here, denotes the polarization vector of the photon with q'^e^ = 0, whereas 
A, V stand for two Lorentz-invariant amplitudes which occur in the general de- 
composition of the tensors 
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I^^ = J dxe^'^^+^'^y (O I TV,^^{x)IUM I K+{p)) , I = V,A . (3.4) 

The form factor A {V) is related to the matrix element of the axial (vector) 
current in ( p.4| ). In Appendix ^ we display the general decomposition of A'^'^, 
y/iu ^Qj, 7^ and provide also the link with the notation used by the PDG [10 
and in O, HI. 



The term proportional to in (|3.2|) does not contain unknown quantities - 
it is determined by the amplitude of the non-radiative decay — > This 
part of the amplitude is usually referred to as " inner Bremsstrahlung (IB) con- 
tribution" , whereas the term proportional to H'^'^ is called " structure-dependent 
(SD) part". 

The form factors are analytic functions in the complex W"^ plane cut along 
the positive real axis. The cut starts at = {Mk + 2M^f for A [at = 
{Mx + Mtt)^ for V]. In our phase convention, A and V are real in the physical 
region of Ki2'y decays, 

<W^ < M|. (3.5) 

The kinematics of (spin-averaged) Ki2j decays needs two variables, for which 
we choose the conventional quantities 

X = 2pq/Ml , y = 2ppi/Ml . (3.6) 

In the K rest frame, the variable x (y) is proportional to the photon (charged 
lepton) energy, 

x = 2E^/Mk , y = 2Ei/MK , (3.7) 
and the angle between the photon and the charged lepton is related to x and 
y by 

{l-y/2 + A/2){l-y/2-A/2) r^—- 
" = l-y/2 + A/2cose,, ' ^ = V^^ " 4n • (3.8) 

In terms of these quantities, one has 

W' = MUI-x) ; iq' = 0) . (3.9) 

We write the physical region for x and y as 

2^1 <y<l + ri 

l-^{y + A)<x<l-^{y-A) (3.10) 



or, equivalent ly, as 



< X < 1 — ri 

x + j-^< y <l + ri (3.11) 
(1 - x) 



where 



ri = m]/M\ 



K 



1.1 ■ 10-6(/ = e) 
4.6-10-2(/ = /i) 



(3.12) 



3.2 Decay rates 

The partial decay rate is 

dV 



J2 \T\'^dLiPs{p;Pi,Pu,q) 



(3.13) 



spins 



The Dahtz plot density 



d^^ _ Mk ^ 1^,2 
dxdy 2567r3 ^ ' ' 



(3.14) 



spins 



is a Lorentz-invariant function which contains V and A in the following form [|T3[1 



p(a;,i/) 
PiB(a;,i/) 
PsD(a;,y) 
Ant (a;,?/) 



Pib(x, ?/) + psD(a;, y) + PiNT(a;, y) 
^iB/iB(a;,z/) 

As^Ml [iy + A) VsD+ (x, y) + iy- Affs^- (x, y)] (3.15) 
A,^^Mk [{V + A)f,^^+{x,y) + {y- A)f,^^-{x,y)] , 



where 

fm{x,y) 

fsB+ix,y) 
/sD- (a;, y) 

fmT+(,X,y) 

/iNT- (a^; Z/) 

and 



1-y + ri 



+ 2(1 -x)(l -n) 



x^(x + y — 1 — r/) 
[x + y - 1 - ri][{x + y - - x) - n] 
[1-y + n] [(l-x)(l-y) + n] 
1 - ?/ + r; 



2xr;(l - ri) 
X -\- y — 1 — ri 

(3.16) 



x{x + y — 1 — ri) 

l-y + ri 
x{x + y — 1 — ri) 



A, 

As: 



[(1 - x){l - X - y) + ri] 
x'^ — (1 — x){l — X — y) — ri 



Ar, 



y2 IT/ |2 



A. 



/Fx 



(3.17) 



4r/ 
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Table 2: The quantities Xj, Nj. 80=*= and INT^ are evaluated with full phase 
space, IB with restricted kinematics ( p.20|) . 





SD+ 


SD- 


INT+ 


INT- 


IB 






1.67-10-2 


1.67-10-2 


-8.22 - 10-® 


3.67-10-6 


3.58 - 10-6 




Xi 


1.18 - 10-2 


1.18 - 10-2 


-1.78- 10-3 


1.23 - 10-2 


3.68- 10-2 




Ni 


2 


2 


1 


1 








For later convenience, we note that 

^SB = ^ ^ J ^)\{K lu,) . (3.18) 
87r ri{l - riY \ Fk J 

The subscripts IB, SD and INT stand respectively for the contribution from inner 
Bremsstrahlung, from the structure-dependent part, and from the interference 
term between the IB and the SD parts in the amplitude. 

To get a feeling for the magnitude of the various contributions IB, SD^ and 
INT^ to the decay rate, we consider the integrated rates 

Tj= [ dxdypi{x,y) ; / = SD^, INT^, IB , (3.19) 

jRi 

where psD = Psd+ + Psd- etc. For the region Rj we take the full phase space for 
/ ^ IB, and 

i?iB = 214.5 MeV/c < pi < 231.5 MeV/c (3.20) 

for the Bremsstrahlung contribution. Here pi stands for the modulus of the lepton 
three-momentum in the kaon rest system Q We consider constant form factors 
V, A and write for the rates and for the corresponding branching ratios 

r, = A,,{MK{V±A)f'Xi 
BRj = Ti/T,^, = N{MKiV±A)}'''Xj (3.21) 

with 

N = Asu/T,^, = 8.348 - IO-2. (3.22) 

The values for Nj and Xj are listed in Table ^. 

To estimate Fj and BR/, we note that the form factors V, A are of order 

Mk{V + A) ^ -10-^ , M/^(V- A) ~ -4- 10-2 . (3.23) 

This cut has been used in for because this kinematical region is free from K^^^ 

background. We apply it here for iUustration also to the electron mode Ke2-y 



9 



From this and from the entries in the table one concludes that for the above 
regions i?/, the interference terms INT^ are negligible in -ft'e27; whereas they are 
important in K^2'y Furthermore, IB is negligible for -ft'e27, because it is helicity- 
suppressed, as can be seen from the factor in Ajb. This term dominates how- 
ever in K^j,2-i- 

3.3 Determination of A{W'^) and V{W'^) 

The decay rate contains two real functions 

F^iW"^) = ViW"^) ± (3.24) 

as the only unknowns. The density distributions fi^, . . . , /int± have very different 
Dalitz plots [0, Therefore, in principle, one can determine the strength of 



each term by choosing a suitable kinematical region of observation. To pin down 
F"^ , it would be sufficient to measure at each photon energy the interference term 
INT^. This has not yet been achieved, either because the contribution of INT^ 
is too small (in -^"627), or because too few events have been collected (in -K'^27)- 
On the other hand, from a measurement of SD^ alone one can determine A, V 
only up to a fourfold ambiguity: 

SD± ^ {{V, A); -{V, A); {A, V); -{A, V)} . (3.25) 

In terms of the ratio 

7K = A/V (3.26) 

this ambiguity amounts to 

SD^ ^ {jk-Mik} . (3.27) 

Therefore, in order to pin down the amplitudes A and V uniquely, one must 
measure the interference terms INT^ as well. 

3.4 Previous experiments 



The PDG uses data from two experiments |[TT|, [T^ , both of which have been 



sensitive mainly to the SD"*" term in ( |3.16| ). In [|T^, 56 events with > 100 



MeV, Ee+ > 236 MeV and 9e+.y > 120° have been identified, whereas the later 



experiment has collected 51 events with > 48 MeV, Ee+ > 235 MeV, 
and 6'e+^ > 140°. In these kinematical regions, background from —>■ e'^UeTi^ is 
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absent because i?™'"'(ii'e3) 

isQ in 



228 MeV. The combined result of both experiments 



r(SD+)/r(ir^2) = (2.4 ± 0.36) ■ 10"^ (3.28) 

For SD^, the bound 

r(SD-)/r,„,,, < 1.6 ■ 10-^ (3.29) 

has been obtained from a sample of electrons with energies 220 MeV < < 230 
MeV jni. Using ( p:2Tl , p:22|) , the result (g) leads to 

Mk\V + A\= 0.105 ± 0.008 . (3.30) 



The bound ( p.29|) on the other hand implies fTT 



V - A \ / \ V + A\<^ll 



(3.31) 



from where it is concluded |TT|] that is outside the range —1.86 to —0.54, 



^[-1.86,-0.54] . 



(3.32) 



As we already mentioned, the interference terms INT^ in — > evel are small 
and can hardly ever be measured. As a result of this, the amplitudes A, V and 
the ratio 7i<- determined from -ft'e27 are subject to the ambiguities (|3.25|) , ( |3.27| ). 



Here, the interference terms INT^ are non-negligible in appropriate regions 
of phase space, see |T3| , |T^. Therefore, this decay allows one in principle to pin 
down V and A. The PDG uses data from two experiments ||12|, |1^. In |jl2[, the 
momentum spectrum of the muon was measured in the region (|3.20| ). In total 
2 ± 3.44 SD+ events have been found with 216 MeV/c < < 230 MeV/c and 

> 100 MeV, which leads to 



Mk\V + A\< 0.16 



(3.33) 



In order to identify the effect of the SD terms, the region 120 MeV/c < < 150 
MeV/c was searched. Here, the background from ii"^3 decays was very serious. 



The authors found 142 Kfj,^,^ candidates and conclude that 



1.77 < Mk{V - A) <0.21. 



(3.34) 



The result ( |3.33|) is consistent with (|3.30|) , and the bound ( p.34|) is worse than 
the result ( ^.31| ) obtained from Ke2-y The branching ratios which follow from 
(|3.33| , |3.34| ) are displayed in Table ^ where we also show the Ke2-y results |jT5|, ^ 
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The entry SD +INT for is based on additional constraints from Kf,2^ 

Hn all four experiments ^ |l6j discussed here and below, the form factors A and V 

have been treated as constants. 
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The Ke2'y data are from 



Table 3: Measured branching ratios V[K lui'y)/Vt^ 
[|T5| , p]T| , the K^2'y data from [|l^, |T^. The last column corresponds to the cut 





SD+ 


SD" 


INT+ 


SD" + INT" 


total 




(1.52 ±0.23) • 10-^ 


< 1.6- 10--* 








K^2-i 


< 3 • 10-^ 




< 2.7 • 10-'^ 
(modulus) 


< 2.6 • 10-'' 
(modulus) 


(3.02 ±0.10) • 10-^ 



3.5 Theory 

The amplitudes y4(PF^) and l^(iy^) have been worked out in the framework of 
various approaches such as current algebra, PCAC, resonance exchange, disper- 
sion relations, .... For a rather detailed review together with an extensive list of 
references up to 1976, see |T^. Here, we concentrate on the predictions of V, A in 
the framework of CHPT. 



A) Chiral expansion to one loop 



At leading order in the low-energy expansion, one has Q 

A = 1/ = (3.35) 

and Fk replaced by F in Eq. (|3.2|) . The rate is therefore entirely given by the IB 
contribution at leading order. The loop effects manifest themselves only in the 
replacement F Fk in Eq. ( p.2|) . The local terms of order give 

V = -^1 

87r2F 

= 32n'iLl + L[,), (3.36) 

where Lg and L\q are the renormalized low-energy couplings evaluated at the 
scale /i (the combination Lg + L[q is scale-independent). The vector form factor 
stems from the Wess-Zumino term 0. 

Remarks : 

^ The relevant Feynman diagrams are displayed and discussed in the next section in con- 
nection with the decays l^vl'^l'^ where the photon is virtual. 
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Table 4: Chiral prediction at order for the branching ratios r{K lvi'-))/T- 
The cut used in the last column is given in Eq. ( |3.2CI| ). 





SD+ 


SD- 


1NT+ 


INT- 


total 




1.30 ■ 10-5 


1.95 ■ 10-6 


6.64-10-1° 


-1.15 ■ 10-s 


2.34 ■ 10-6 




9.24 ■ 10-6 


1.38 ■ 10-6 


1.44- 10-5 


-3.83-10-5 


3.08 ■ 10-3 



1. At this order in the low-energy expansion, the form factors A, V do not exhibit 
any W"^ dependence. A non-trivial W"^ dependence only occurs at the next 
order in the energy expansion (two- loop effect, see the discussion below). Note 



that the available analyses of experimental data of — >^ lui'j decays [1^, |TT 



12, 16] use constant form factors throughout. 



Once the combination Lg + Liq has been pinned down from other processes, 
Eq. (p.36| ) allows one to evaluate A, V unambiguously at this order in the 
low-energy expansion. Using Lg + Liq = 1.4 ■ 10-^ and F = F^, one has 

Mk{A + V) = -0.097 
MKiV-A) = -0.037 

7x = 0.45 . (3.37) 



The result for the combination [A + V) agrees with ( |3.30| ) within the errors, 
while 7x is consistent with ( p.32|) . 

3. At this order in the low-energy expansion, the form factors A, V are identical 
to the ones in radiative pion decays 7r+ —>■ I'^ui'-f 0, as a result of which one 
has 

ln = lK + OiM) , (3.38) 



We display in Table ^ the branching ratios BR/ (|3.21| ) which follow from the 



prediction ( |3.37|) . These predictions satisfy of course the inequalities found from 



experimental data (see Table H). 

B) W'^ dependence of the form factors 

The chiral prediction gives constant form factors at order p'^. Terms of order 
p6 have not yet been calculated. They would, however, generate a non-trivial 
dependence both in V and A. In order to estimate the magnitude of these 
corrections, we consider one class of contributions: terms which are generated 
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Figure 1: The rate dP{x)/dx in ( |3.40|) , evaluated with the form factors ( |3.39| ) 
and N,,, = 9 • 10^. The sohd hue corresponds to Mk* = 890 MeV, Mk, = 1.3 
GeV. The dashed hne is evaluated with M^* = 890 MeV, M^^ = oo, and the 
dotted line corresponds to M^* = = oc. The total number of events is also 
indicated in each case. 



by vector and axial vector resonance exchange with strangeness |17, 18 



V(W^) = ^1 ^ , A(W^) = 4 ^ , (3.39) 

where A are given in ( |3.36|) . We now examine the effect of the denominators 
in (|3.39|) in the region y > 0.95, x > 0.2 which has been explored in — > e~^Ue^ 



np. We put rUe = and evaluate the rate 

dP{x) N,,, 



, -p , PsD+ix,y)dy , (3.40) 

dx 1 tot Jy=0.95 

where N^^^ denotes the total number of decays considered, and T'^J: = 1.24 ■ 
10-8 s. 

The function is displayed in Fig. ^ for three different values of Mk* and 
Mki, with Ntot = 9 ■ 10^. The total number of events 

Np= dP{x) (3.41) 

Jx=Q.2 

is also indicated in each case. The difference between the dashed and the dotted 
lines shows that the nearby singularity in the anomaly form factor influences the 
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decay rate substantially at low photon energies. The effect disappears at x ^ 1, 
where W"^ = M]^{1 — x) — > 0. To minimize the effect of resonance exchange, the 
large-x region should thus be considered. The low-x region, on the other hand, 
may be used to explore the dependence of V and of A. For a rather exhaustive 
discussion of the relevance of this dependence for the analysis of Ki2'y decays 



we refer the reader to Ref. \17 



3.6 Outlook 

Previous experiments have used various cuts in phase space in order (i) to identify 
the individual contributions IB, SD^, INT^ as far as possible, and (ii) to reduce 
the background from Ki^ decays. This background has in fact forced such severe 
cuts that only the upper end of the lepton spectrum remained. 

The experimental possibilities to reduce the background from Ki^, decays are 
presumably more favourable with today's techniques. Furthermore, the annual 
yield of 9 ■ 10^K~^ decays at e.g. DAFNE is more than two orders of magnitude 



higher than the samples which were available until now [|TT| |T2|, |T5|, [T^. This 
allows for a big improvement in the determination of the amplitudes A and V, in 
particular in decays. It would be very interesting to pin down the combina- 
tion L9 + L10 of the low-energy constants which occurs in the chiral representation 
of the amplitude A, and to investigate the dependence of the form factors. 



4 The decays K^, l^uil'+l'- 

Here we consider decays where the photon turns into a lepton-antilepton pair. 



- 






(4.1) 


K+ - 






(4.2) 


7r+ - 






(4.3) 


- 






(4.4) 


K+ - 






(4.5) 


7r+ - 






(4.6) 



4.1 Matrix elements 

We start with the processes ( [4.1D and (|4^ ), 



(/,/') = (e,/x)or (/i,e). (4.7) 

The matrix element is 

T = -iGFeV:rep [FkL' - Wl,} (4.8) 



where 



2pq — Zpiq + q 



= iVitP^^'^q^pp- A^^qWgP^" -WPq^") 

-yl2(gV^ - q^q^) - Ai{qWqP - q^WP)W^' , (4.9) 



with 



.,2 \n J +-^3- (4.10) 



2Fk F^jq' ) - 1 
The form factors ^^(g^Vr^)^ y^(^q2^^^2^ 

are the ones defined in Appendix B, 

and Fy{q^) is the electromagnetic form factor of the . Finally the quantity 
stands for 

= ^u{p2)rv{p^) , (4.11) 

and the four-momenta are 

q = Pi+P2, W=pi+pi, = p-q, (4.12) 



such that q^I^ = 0. The first term in Eq. (|4.8| ) contains the part where the 
off-shell photon radiates off the final-state lepton and those corrections that are 
reabsorbed in the definition of Fx- This can always be done independently of the 
model used, because of gauge invariance [0. 

In order to obtain the matrix element for ([4.4|) and ([4 .51), 



K+{p) ^ l^{piMp.)l^{pi)l-{p2) , (4.13) 



one identifies mi and mj in ( [4.8|) and subtracts the contribution obtained from 
interchanging pi ^ pi : 

(puPi) (PhPl) 

q P1+P2 

W -> p-q = p^+pi . (4.14) 

For the decays of a pion, the same formulas with the following replacements 
apply: 

Mk^M^, Vus Vud , 

Fk-F^, F^iq') ^ F^iq') . (4.15) 
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4.2 Decay distributions 

The decay width is given by 



= nnj \o ^8 X! \T^\dLiPs{p;Pi,Pu,Pi,P2) (4.16) 



2Mi^(27r)8 



spins 



and the total rate is the integral over this for the case / 7^ /'. For the case / = /' 
the integral has to be divided by a factor of 2 for two identical particles in the 
final state. 

We first consider the case where I 7^ V and introduce the dimensionless vari- 
ables 

_ 2m _ 2pip _ 

^ ~ Ml- ' Ml ' Ml ' 

Then one obtains, after integrating over pi and p2 at fixed 



spins 



, „ 1 ( 2r,'| 4r; 



T" = Af^2{FKL''-ff''%}. (4.18) 

The quantity |— Z!spms ^^^^} can be found in Ref. Jl^. The expression or a 
Monte Carlo program containing this trace in FORTRAN can be obtained from 
the authors. This result allows one to evaluate, e.g. the distribution dV/dz of 
produced I'^l'^ pairs, rather easily. The kinematically allowed region is 

Ar[< z < 1 + n - 
2^ < X <l + z-ri (4.19) 
A- B < y <A + B, 

with 

[2 - x){l + z + ri - x) 



A 



2(l + ;z-a;) 



_ {l + z-x-ri)^x'^ -Az 

2(l + ;z-a;) ■ ^ ' ' 

The case / = /' is slightly more elaborate since the integration over part of the 
phase space cannot be done analytically as before. The expression for Y.snins I^P) 
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Figure 2: The two classes of contributions to the decays — > l^vil'^l' . 

together with the Monte Carlo program to do the phase-space integrals, is avail- 
able from the authors. 

All formulas in this section are also valid with the replacements of Eq. (|4.15| ) 
for the pion-decay case. 



4.3 Theory 

The form factors Aj, Vi and Fy have been discussed in all kinds of models, 



Vector Meson Dominance, hard meson, etc. (see Ref. for a discussion). Here, 
we will work rigorously within the framework of CHPT. 

There are two classes of contributions to these decays. They are depicted in 
Fig. 1^. The first class of diagrams is the well understood Bremsstrahlung off the 
final-state lepton. This contribution is entirely Q contained in the first term in Eq. 
( [4.8|) . The second class of diagrams, which also contains the radiation off the kaon 
line, contributes to both terms in Eq. ( |4.8| ). The one- loop diagrams contributing 
to this class are shown in Fig. 

To leading order we have 

1^1 = 

= A2 = A3 = . (4.21) 

We also have Fy = 1 and F^ is replaced by F. The rate here is entirely given 
by the inner Bremsstrahlung contribution. 

^The momentum dependence of the W propagator is neglected. 
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Figure 3: Loop diagrams for — >■ The virtual photon must be appended 

on all charged lines and on all vertices. 

The effects of the next-to-leading order CHPT corrections manifest themselves 
in the replacement F —>■ Fx and in non-zero values for several of the form factors 



1 



87r2F 

Ai = -j{Ll + L[,) 

A. ^ JJdfim^ (4.22) 
A3 = 



These results obey the current algebra relation of Ref. The function Fy{q^) 
does, however, deviate somewhat from the linear parametrization that is often 
used. The function H{t) is defined in Appendix 0. It can be easily seen from ( ^.8| ), 
L9|) and ([4.23|) that the one-loop contributions to the decay matrix elements 



vanish for — * and only the contributions from the Wess-Zumino term (to 
the form factor Vi) and from the Lagrangian remain. 

The fact that the form factors at next-to-leading order could be written in 
terms of the kaon electromagnetic form factor in a simple way is no longer true 
at the level. The Lagrangian at order contains a term of the form 

ti {D^FrU^D^FRp.u} (4.23) 

that contributes to A2 and ^3, but not to the kaon electromagnetic form factor. 
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Table 5: Theoretical values for the branching ratios for the decay K 
^'^h'^e~^e~ for various cuts. 





tree level 


form factors as given by CHPT 


full phase space 


2.49 ■ 10-5 


2.49 ■ 10-5 


z < 10-3 


2.07-10-5 


2.07- 10-5 


z > 10-3 


4.12 ■ 10-6 


4.20 ■ 10-6 


z > (20 McV/Mk)'^ 


3.15-10-6 


3.23 ■ 10-6 


z > (140 McV/Mk)'^ 


4.98 ■ 10-8 


8.51 ■ 10-8 


X > 40 McV/Mk 


1.58 ■ 10-5 


1.58 • 10-5 



Essentially the same result holds for the decay 7r+ /+z/;e+e- after doing 
the replacements ( [4. 151 ) in Eqs. ( |4.8| ), ( [4. 91 ) and ( [4.231 ). The pion electromagnetic 
form factor to one loop is given by pT| 

F^{q^) = l + 2H^^{q^) + HKK{q^) • (4.24) 

Our result for the vr"*" — ^ l^vil^l^ form factors agrees in the limit of two flavours 
with the result of Ref . . 



4.4 Numerical results 

Using the formulas of the previous subsections, we have calculated the rates for 
a few cuts, including those given in the literature. The values for the masses used 
are those of and vr"*". For Lg and Lio we used the values given in Section ^. 

For the case of unequal leptons, the results are given in Table ^ for the 
decay —>■ ii^v^e'^e~. These include the cuts used in Refs. and [^ , 



X > 40 MeV/Mx and z > (140 MeV/MK)^ respectively. It can be seen that 
for this decay most of the branching ratio is generated at very low electron- 
positron invariant masses. As can be seen from the result for the cuts used in 
Ref. [^, the effect of the structure-dependent terms is most visible at high in- 



variant electron-positron mass. Our calculation, including the effect of the form 
factors, agrees well with their data. We disagree, however, with the numerical 
result obtained in Ref. by about an order of magnitude. 

For the decay — > e"'"z/e/i "*"//-, we obtain for the tree level or IB contribution 
a branching ratio 

BRib{K+ e+i^e/^V) = 3.06 ■ 10-^^ (4.25) 
and, including the form factors, 

BRtot.i{K+ -> e+z/e/iV~) = 1-12 ■ 10"^ (4.26) 
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Table 6: Theoretical values for the branching ratios for the decay 7r+ ji^u^e^e 
for various cuts. The effect of the form factors is not visible to the accuracy quoted. 





tree level 


full phase space 


3.3-10-^ 


z > 10-3 


6.2 ■ 10-8 


Z > 10-2 


2.8 ■ 10-9 



Here the structure-dependent terms are the leading contribution since the inner 
Bremsstrahlung contribution is helicity-suppressed, as can be seen from the factor 
m in L^. 

The decay vr^ n^u^e^e" is entirely given by the inner Bremsstrahlung 
contribution. The effects of the form factors are very small. The IB amplitude is 
not helicity-suppressed here since /M^ pa 1 . The effects of the form factors are 
partly suppressed since only very small e~^e~ masses are allowed by phase space. 
The rates for full phase space and various cuts on z, defined as in Eq. (|4.29|) , but 
with instead of M^, are given in Table ^. 

For the decays with identical leptons we obtain for the muon case a branching 
ratio of 

BRtoUK+ ^ /i+z/^/iV") = 1-35 ■ 10-^ (4.27) 

for the full phase space, including the effects of the form factors. The inner 
Bremsstrahlung or the tree level branching ratio for this decay is 

BRib{K+ = 3.79 • 10-9. (4.28) 

For the decay with two positrons and one electron, the integration over the full 
phase space for the tree-level results is very sensitive to the behaviour for small 
pair masses. We have given the tree level and the full prediction including form- 
factor effects in Table |^. The cuts are always on both invariant masses : 

zi = {pi+P2?/Ml. (4.29) 



For the decay 7r+ e+i/eC^e-, the integration over phase space is again very 
sensitive to the region for small invariant masses of the pairs. Tree level and full 
results are given for several cuts in Table ||; z and zi are defined as in Eq. ( |4.29| ) 
with Mk replaced by M^. The cuts used for the last row are similar to those of 
Ref. |2^. We require a minimum energy of 15 MeV for the three charged particles 



and an opening angle 6i 2 for each electron-positron pair with cos^i 2 < 0.95. 
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Table 7: Theoretical values for the branching ratios for the decay e^u^e^e 
for various cuts. 





tree level 


form factors as given by CHPT 


full phase space 


^ 4 ■ 10-9 


1.8 ■ 10-^ 


z, zi > 10"^ 


3.0 ■ 10-10 


1.22 ■ 10-'^ 


z, zi > (50 MeV/MKf 


5.2 ■ 10-11 


8.88 ■ 10-s 


z, zi > (140 MeV/MKf 


2.1 • 10-12 


3.39 ■ 10-^ 



Table 8: Theoretical values for the branching ratios for the decay vr''" e^Uee^e 
for various cuts. 





tree level 


form factors as given by CHPT 


z, zi > 10-3 


2.4- 10-9 


3.0-10-9 


Z, Zi > 10-2 


4.2-10-10 


8.7-10-10 


Cuts similar to 23] 


2.1 ■ 10-10 


5.7-10-10 



4.5 Experimental status 

For kaon decays only those with an electron positron pair in the final state, decays 
( [4.2| ) and ([4.4|), have been observed. 

Both were measured in the same experiment [^. The decay ^^u^e^e^ 
was measured with a branching ratio of (1.23 ± 0.32) - 10-^ with a lower cut on 
the electron-positron invariant mass of 140 MeV. The measurement is compatible 
with our calculation including the form-factor effects for the relevant region of 
phase space. This measurement was then extrapolated p2| using the result of pO 



to the full phase space. Since we disagree with that calculation, we also disagree 
with the extrapolation. 

In the same experiment, four events of the type e'^u^e'^e' were 

observed where both electron-positron pair invariant masses were above 140 
MeV. This corresponds to a branching ratio for this region of phase space of 
{2.8tiA) ' IQ-^. This result is compatible within errors with our calculation, see 
Table 0. The matrix element of Ref. |2^ was again used for the extrapolation to 
full phase space ||2^. Apart from our numerical disagreement, the calculation of 
Ref. ||20[ was for the case of non-identical leptons and cannot be applied here. 

For the decay —>■ fi'^u^^fi'^ fi~ an upper limit of 4.1 - 10-'' exists [^]. This 
upper limit is compatible with our theoretical result, Eq. ( [1.27| ). 



The decay e'^u^ii^n has not been looked for so far but should be 

within the capabilities of a $ factory, given the branching ratio predicted in the 
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previous subsection. This decay proceeds almost entirely through the structure- 
dependent terms and, as such, is a good test of our calculation. 

The decay 7r+ fi~^i'^e~^e~ has not yet been looked for. It is almost entirely 
given by the inner Bremsstrahlung contribution so it is not very interesting as a 
test of the form factors. The decay tt~^ —>■ e^Vf.e^e~ has been measured in 
We do not make a direct comparison with our result, because the acceptance of 



the apparatus is very difficult to implement, see Ref. ||26|. We do agree with the 
form factors as determined there. 



4.6 Outlook 

The decays discussed in this section, l^uil'^l'^ , are complementary to the 

decays — >• I'^i'i'j. As was the case for the analogous decay 7r+ e^u^e^e~ 
[ P^ , it may be possible to explore phase space more easily with this process than 
with l^vn to resolve ambiguities in the form factors. 

As can be seen from our predictions. Tables H and |^, all the decays considered 
in this section should be observable at a $ factory like DAFNE or in fixed-target 
experiments at future kaon factories. Large improvements in statistics are possible 
since less severe cuts than those used in the past experiments should be possible. 
In the decays with a fi^fi^ pair and the decay e^u^e^e^ , the effects of the 

form factors are already large in the total rates. In the decay yU+z/^e^e~ 
most of the total rate is for a small invariant mass of the pair and is given by the 
inner Bremsstrahlung contribution. There are, however, regions of phase space 
where the form-factor effects are large and where branching ratios are still large 
enough that sufficient statistics can be accumulated. 



5 Radiative Kis decays 

The decay channels considered in this section are 

K\p) n-{p'np,MpMQ) 
and the charge conjugate modes. We only consider real photons (g^ = 0). 
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5.1 Kinematics and invariant amplitudes 

The matrix element for has the general structure 



T = ^eK>"(g)*|(y+-A-MP.)7'^(l-75)^(p/) (5.1) 

2piq J ^ 

The diagram of Fig. corresponding to the first part of Eq. ( |5.1|) , includes 
Bremsstrahlung off the K~^. The lepton Bremsstrahlung diagram of Fig. ^ is 
represented by the second part of Eq. (|0|). The hadronic tensors V^^.A^^ are 
defined as 

/+ = I J d'xe'^^TT^p') I T{V^-ix)lt-''{0)} I K+ip)) , I = V,A; (5.2) 
is the matrix element 

F;^ = (Ap') I K^'^'(O) I K^p)) . (5.3) 
The tensors and satisfy the Ward identities (cf. Appendix C) 

q'v;, = K (5.4) 

q^A-^^, = 0, 

leading in turn to 

gM+ = , (5.5) 

as is required by gauge invariance. 

For -R'/s^, one obtains the corresponding amplitudes and hadronic tensors by 
making the replacements 

'^X ~^ 1 ^ A% 5 (5.6) 

To make the infrared behaviour transparent, it is convenient to separate the 
tensors V^, V^^ into two parts: 

= V;t + ^F+, (5.7) 
^ pq 
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Due to Low's theorem |]T9|, the amphtudes V^^'° are finite for g — > 0. The axial 
amphtudes ^^j^^ are automatically infrared- finite. The Ward identity ( ^.41) implies 
that the vector amplitudes V^'^ are transverse: 



(5i 



For on-shell photons, Lorentz and parity invariance together with gauge in- 
variance allow the general decomposition (dropping the superscripts +, and 
terms that vanish upon contraction with the photon polarization vector) 



V, 



Vi gf,u 



qW ^ 



(5.9) 



pq 



+V3 I p',w. - ^w,w, ]+V,{ p'X 



qW 



pq 
qW 



V = le^^pMiP'^q" + A^qW) + ie^xpaV'\'W''{AzW, + A^p'J) 
= C^p', + C2{p-p')u 

W = Pl+Pu ■ 

For the decomposition of the axial tensor amplitude, use was made of Schouten's 
identity 

VxSf^vpa = , (5.10) 

where the sum extends over all cyclic permutations of the five indices and v is 
an arbitrary four-vector. With the decomposition ( |5.7D we can write the matrix 
element for i^T^g^ in ( |5.1| ) in a form analogous to Eqs. (|3?^ , |37^ ) for Ki2^: 



T 



V2 

+ F+n(p,)7'^(l-75 



(5.11) 



pq 



^piq 



The four invariant vector amplitudes , . . . , V4 and the four axial amplitudes 
Ai, . . . , ^4 are functions of three scalar variables. A convenient choice for these 
variables is 

E^ = pq/MK, E^=pp'/Mk, W = VW^ , (5.12) 

where W is the invariant mass of the lepton pair. The amplitudes Ci,C2 can 
be expressed in terms of the Ki^ form factors and depend only on the variable 
(p-p')2 = MI + MI-2MkE^. For the full kinematics oiKi^^ two more variables 
are needed, e.g. 

El = ppi/Mk, X = piq/Ml . (5.13) 
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Figure 4: Diagrammatic representation of the amplitude. 

The variable x is related to the angle between the photon and the charged 
lepton in the K rest frame: 

xMl = E^{Ei - ^Ef -mf cos Oi^) . (5.14) 

T invariance implies that the vector amplitudes ^i, . . . , V4, the axial amplitudes 
y4i, . . . , and the Ki^ form factors Ci, C2 are (separately) relatively real in the 
physical region. We choose the standard phase convention in which all amplitudes 
are real. 

For Oi-y -—>■ (collinear lepton and photon), there is a lepton- mass singularity 
in (|5.1|) , which is numerically relevant for / = e The region of small E.y, 61^ is 
dominated by the K13 matrix elements. The new theoretical information of Ki^^ 
decays resides in the infrared-finite tensor amplitudes V^u and A^i,. The relative 
importance of these contributions can be enhanced by cutting away the region of 

low Ey, 9ly. 

5.2 CHPT to 0(/) 

Prior to CHPT, the most detailed calculations of Ki^^ amplitudes were performed 
by Fearing, Fischbach and Smith using current-algebra techniques. 

To lowest order in the chiral expansion [0(p^)], the invariant amplitudes for 
the two charge modes are given by |@] 
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Figure 5: Loop diagrams (without tadpoles) for Ki^,. For -ft'ia-y, the photon must 
be appended on all charged lines and on all vertices. 



= ^ (5.15) 

V2pq 

C+ = 2C+ = V2 



= -1 

= ;^ (5.16) 



p'q 



All other invariant amplitudes vanish to O(p^). 

At next-to-leading order in CHPT, there are in general three types of con- 
tributions 0: anomaly, local contributions due to £4, and loop amplitudes. The 
chiral anomaly encoded in the WZW functional yields the axial amplitudes 



iV2 



= I + Wr + ^,_ ^.^ e^.^MYW^j (5.17) 
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Table 9: Coefficients for the K^^ loop amplitudes corresponding to the diagrams 
J = 1, 2, 3 in Fig. |[ All coefficients c[ must be divided by 6v^F^. 



/ 


Mi 


mi 


c{ 






1 


Mk 


M^ 


1 


-2 


-M|. - 2M2 


2 


Mk 


Mr, 


3 


-6 


-Ml - 2Ml 


3 


M^ 


Mk 





-6 


-ml 



The local parts of these amplitudes can be read off the Lagrangian ( p. 181 ). The 
non-local kaon pole term in ^4+^ is due to the anomalous K^K'n^'-y vertex con- 
tained in ( p.l9| ) with the outgoing turning into a W'^ . Its contribution to 
the matrix element is helicity-suppressed p8| , |^ and can be neglected for / = e. 
However, for K'^-^^ it is a normal-size O(p^) contribution to the axial amplitude. 

The loop diagrams for Ki^^ are shown in Fig. |^. We first write the K'^ matrix 
element in terms of three functions /{*", f2-, ft which will also appear in the in- 
variant amplitudes V^^ . Including the contributions from the low-energy constants 
L5, Lg in £4, the matrix element F+ is given by 



ct 
ct 

ft it) 
ftit) 
ft it) 

A/ 



ft it) 

\iMl-Ml-t)ftit) + ftit) 



V2 + 
4L9 



AT 

-^^t + 2j:ici-ci)Biit) 




2Bi{t) 



{t + A,)Aj J,(t) 
2t 



^ + 4 E [i^i + + A/) - H) ^iJiit) 
M^MlM,, 




2567r2 



In- 



Mi 



m 



I 5 



t 



ip-p 



i\2 



Lg is a scale-independent coupling constant and we have traded the tadpole con- 
tribution together with L5 for Fk/F^ in ft it) [cf. Eq. ([2T^ )]. The sum over I 



corresponds to the three one- loop diagrams of Fig. |], with the coefficients c(, C2, Cg 
displayed in Table ^ We use the Gell-Mann-Okubo mass formula throughout to 
express M^ in terms of M^^M^. The functions Jj{t) and B^{t) can be found in 
Appendix 0. 
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The standard Ki^ form factors f+(t), f-(t) are 
1 

1 

71 



hit) 
f~{t) 



ftit) 

{Ml-Ml-t)ft{t)+2ft{t)-ft{t) 



(5.19) 



It remains to calculate the infrared-finite tensor amplitude V^. The invariant 
amplitudes V^^ can be expressed in terms of the previously defined functions f^' 
and of additional amplitudes Ii, 12,13. Diagrammatically, the latter amplitudes 
arise from those diagrams in Fig. |^ where the photon is not appended on the 
incoming (non-Bremsstrahlung diagrams). The final expressions are 



h+p'wjnw;) + f^{w;) 



pq 
1 



V^^ = h + - p'WftiW') + f+{W') - p'WJ+{W') - f+{W') (5.20) 



pq 



v; 



pq 

W + q = p — p' . 



The amplitudes /i, /2, h in Eq. ( ^.20| ) are given by 
+ E ( [iK + A,)(c{ + ci) - 2{cip'W, + ci) 



1=1 



{ci - ci) 



mi 



p'W,AW^-Aj)Jj 



q 



ABi)+p'{W-q)Ll 



^^-^ [p'q{F, - {Wl + A,)G,) + p'W{Bi - Bi)] } 



h 
h 

lL 



qW 

+47j:{ci-c{)\Fj-{W' + Aj)Gj 



v/2F2 qW fr{ 



4Lq 



V2F^ 



+ E 2(c^-c{ 

7=1 I 

M^MlM, 



^1 



2567r2 



In- 



(5.21) 



29 



Table 10: Coefficients for the /^^^^ loop amplitudes corresponding to the diagrams 
/ = 1, 2, 3 in Fig. ^. All coefficients cf must be divided by 'o\f2F'^ . 



I 


Mi 


nil 


c{ 






1 


Mk 







-3 


-'iMl 


2 


Mk 


Mr, 


6 


-3 


Ml + 2Ml 


3 




Mk 


4 


-2 


-2M| + 2Ml 



Gi 
Ji 



W^Bi - W^Bi 



{W^^+/\i)Ji-{W^ + /\i)Ji 



+ 



647r2 



The function C{Wq,W'^, Mfj-mj) is given in Appendix 0. All the invariant 
amplitudes V^, . . . , are real in the physical region. Of course, the same is true 
for the Ki^ form factors , C2 ■ 

The Kf^^ amplitude has a very similar structure. Both the Kf^ matrix ele- 
ment and the infrared-finite vector amplitude V^^^ can be obtained from the 
corresponding quantities F+ and by the following steps: 

• interchange p' and —p ; 

• replace ^ by ^ in /3+; 

• insert the appropriate coefficients cf for Ki^^ listed in Table 

• multiply and by a factor — -\/2. 



5.3 Numerical results 

The total decay rate is given by 
r{K Tclui'-f) - 



2Mk{2'kY 



JdLiPs{p;p',Pi,Pu,q) E 1^1' (5-22) 



spms 



in terms of the amplitude T in (|5.1| ). The square of the matrix element, summed 
over photon and lepton polarizations, is a bilinear form in the invariant ampli- 
tudes Vi, . . . ,Vi, Ai, . . . , A4^,Ci, 6*2. The explicit expression can be found in Ref. 
TM, which also contains a summary of the presently available data sample. 
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Table 11: Branching ratios and expected number of events at DAFNE for Kf 





BR 


^events/yr 


full O(p^) amplitude 
tree level 

O(p^) without loops 


3.0 X 10-^ 
2.8 X 10-"^ 
3.2 X 10"^ 


2.7 X 10^ 
2.5 X 10^ 
2.9 X 10^ 






BR 


#events/yr 


full 0{p'^) amplitude 
tree level 

O(p^) without loops 


2.0 X 10-5 
1.9 X 10-5 

2.1 X 10-5 


1.8 X 105 
1.7 X 105 

1.9 X 105 



Detailed experimental studies of Ki^^ decays will become feasible with the 
high kaon fluxes of various factories, which are either under construction or still 
in the planning stage. The following number of events corresponds to the kaon 
fluxes expected in the proposed detector KLOE for the $ factory DAFNE in 



Frascati. These fluxes are based on a luminosity of 5 x 10'^^ cm^^ s^^ equivalent 
to an annual rate of 9 x 10^ (1.1 x 10^) tagged (-^l) assuming a year of lO'^ 
s. 

For the calculation of decay rates we have used the following cuts: 

> 30 MeV (5.23) 
01^ > 20°. 

The physical values of M^^ and Mk are used in the amplitudes; is calculated 
from the Gell-Mann-Okubo mass formula. The values of the other parameters can 
be found in Sect. |^. For Kf^^, the rates are to be understood as T{Kl n'^Puj). 

The results for K^^ and Kj^^^ are displayed in Tables |1T] and |T2|, respectively. 



For comparison, the tree-level branching ratios and the rates for the amplitudes 
without the loop contributions are also shown. The separation between loop and 
counterterm contributions is of course scale-dependent. This scale dependence is 
absorbed in the scale- invariant constants Lq,Lio defined in Eqs. ( |5.18| ), ( |5.21| ). 



In other words, the entries in Tables 11, 12 for the amplitudes without loops 



correspond to setting all coefficients c[ in Tables equal to zero. 

These numerical results show very clearly that the non-trivial CHPT effects 
of 0(]9^) can be detected at dedicated machines such as DAFNE, in all four chan- 
nels, without any problem of statistics. Of course, the rates are bigger for the 
electronic modes. On the other hand, the relative size of the structure-dependent 
terms is somewhat bigger in the muonic channels (around 8% for the chosen cuts). 
We observe that there is negative interference between the loop and counterterm 
amplitudes. With the chosen distinction between loop and counterterm contribu- 
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Table 12: Branching ratios and expected number of events at DAFNE for Kl 



-"■637 


BR 


#events/yr 


full O(p^) amplitude 
tree level 

0{p^) without loops 


3.8 X 10-3 
3.6 X 10-3 
4.0 X 10-3 


4.2 X 10^ 
4.0 X 10^ 
4.4 X 10^ 



^m37 


BR 


^events/yr 


full 0{p^) amplitude 
tree level 

0{p^) without loops 


5.6 X 10"^ 
5.2 X 10-^ 
5.9 X 10-^ 


6.1 X 10^ 
5.7 X 10^ 
6.5 X 10^ 



tions discussed above, the loops amount to approximately 50% of the rate due to 
the local amplitudes of O(p^) (including the anomalous ones). 

The sensitivity to the counterterm coupling constants Lg, Liq and to the chiral 
anomaly can be expressed as the difference in the number of events between the 
tree level and the O(p^) amplitudes (without loops). In the optimal case of K^^^^, 
this amounts to more than 4 x 10^ events/yr at DAFNE. Most of this difference 
is due to Lg. 

The chiral anomaly is more important for K^^^, but even there it influences 
the total rates rather little. A dedicated study of differential rates is necessary to 
locate the chiral anomaly in Ki^^ amplitudes, if this is at all possible. 

However, taking into account that both Lg and Lio are already known to 
rather good accuracy (see Sect. |]), Ki^^ decays will certainly allow for precise 
tests of CHPT to O(p^). 



6 Conclusions 



In this paper we have calculated the matrix elements for the radiative 
semileptonic kaon decays l~^^n, l^i'il'^l'' and K 7r/i/;7(/, /' = e, /i) 

to next-to-leading order in CHPT. All measured semileptonic K decay 
amplitudes are now known ^ ^ to order in CHPT, except for 



one form factor in K^4 decays ^2 



The Ki2^ decay does not receive any contribution from the one-loop cor- 
rections, except for those that can be absorbed into F^- This decay is then 
completely predicted, including effects of SU{3) breaking, from the knowl- 
edge of the 7 parameter as measured in 7r+ e^Vel- This prediction was 
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successfully confronted with the present experimental data, and the rele- 
vance of more detailed measurements at future kaon facilities was discussed. 

3. For the decays with a lepton-antilepton pair in the final state {K'^{7r'^) — > 
l'^i>il''^l'~), the additional effects of one- loop contributions can be described 

by the electromagnetic form factor. It was shown that this relation breaks 
down at order . We noted that decays with a fi'^fi~ pair are especially 
useful probes of the form factors. All these decays should be measurable in 
the near future with sufficient statistics to measure the form factors. 

4. The last considered decay, Ki^^^, is a useful test of the non-trivial aspects of 
CHPT since the one-loop effects are of the same magnitude as the effects 
of the tree-level O(p^) action. All form factors vary rather significantly over 
the available kinematical region, so that an analysis of experiments with a 
constant form factor will not be able to uncover these effects. We therefore 
urge future experiments to perform the analysis of their data directly with 
the formulas given here. 

5. Overall, the agreement with the presently available data is quite satisfac- 
tory. On the other hand, future kaon facilities like DAFNE will have the 
opportunity to test the predictions of CHPT at next-to-leading order in 
much more detail than is possible with present data. 
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APPENDICES 



A Notation and loop integrals 

The notation for phase space is the one without the factors of 27r. For the decay 
rate of a particle with four-momentum p into n particles with momenta pi, . . . , pn, 
this is 

dLIPs{p;Pl, - ■ ■ ,Pn) = \p-^Pi) n • (A-l) 

\ i=l ) i=l '^Pi 

We use a covariant normalization of one-particle states, 

(p-"b-) = (2vr)W(p"-p) , (A.2) 
together with the spinor normalization 

u{p,r)u{p, s) = 2m6rs ■ (A. 3) 

The kinematical function X{x,y, z) is defined as 

A(x, z) = + + — 2{xy + yz + zx) . (A. 4) 

We take the standard model in the current x current form, i.e. we neglect the 
momentum dependence of the W propagator. The currents used in the text are : 

= g7M^(A4-«A5)g = s-f^u 



2 



Q = diag(2/3, -1/3, -1/3) . (A.5) 

The numerical values used in the programs are the physical masses for the parti- 
cles as given by the Particle Data Group [0]. In addition we have used the values 



for the decay constants derived from the most recent measured charged pion and 



kaon semileptonic decay rates [M, |33 



= 93.2 MeV 
Fk = 113.6 MeV . (A.6) 

We do not need values for the quark masses. For the processes considered in this 
paper we can always use the lowest-order relations to rewrite them in terms of 
the pseudoscalar meson masses (see Section ^). For the KM matrix element \Vu. 



we used the central value, 0.220, of Ref. |jTO|]. The numerical values for the Ll{Mp) 
are those given in Section ^ 
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Whenever we quote a branching ratio for a semileptonic K'^ decay, it stands 
for the branching ratio of the corresponding decay, e.g. 

BR{K^ TT-l+ui) = BR{Kl -n^l^v) . (A.7) 

We use the Condon-Shortley phase conventions throughout. 

Next we define the functions appearing in the loop integrals used in the text. 
First we define those needed for loops with two propagators, mainly in the form 
given in Ref. @]. We consider a loop with two masses, M and m. All needed 
functions can be given in terms of the subtracted scalar integral J(t) = J(t) — J(0), 

, , . r d^p 1 ( \ Q \ 

with t = k"^. The functions used in the text are then : 

, 1 /-i , M2 - tx(l -x)- Ax 

J(t) = -— dx log ^ ' 



167r2 Jo M2 - Ax 

If A, m2 S, m2 (t + ^)2 _ 

' ^ ' — log — r — log — r log 



327r2 \ t ^ A ^ t ^ {t - ^/Xf - A^ 

.Fit) = J{t)-2k, 

M'it) = —{t-2T.}J{t) + ^J{t) + —^-- 
^ ' I2t^ ^ ^ ' 3t^ ^ ' 2887r2 6 

1 f„ M'^m'^ , rn? ^ 



gevr^tl^ + ^^iog^l , 

A2 - 

m = —Jit), 

Kit) = 

Hit) = l^,t + j-,[tM^it)-Lit)], 

A = M'^-m^ , 
S = + , 

A = A(t,M^m2) = (t + A)2 - 4tM=^ . (A.9) 

In the text these are used with subscripts, 

J.^.(t) = Jit) with M = Mi,m = Mj (A.IO) 

and similarly for the other symbols. The subtraction-point-dependent part is 
contained in the constant k 



k = ^, -HJi (A.ll) 
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where fi is the subtraction scale. 

In addition, in Section ^ these functions and symbols appear in a summation 
over loops / with 

Jj[t) = J{t) with M = Mi,m = mi, 

= Mj + m] (A. 12) 

and again similarly for the others. There the combination B2 appears as well : 

B2{t,M\m^) = 52(^,m^M2) (A.13) 

1 , , X(t,M^,m'^)J(t) tE-SMW, M2 
(3S -t) — + — — log — . 



2887r2 ' ' 12t 3847r2A 

The last formula to be defined is the three-propagator-loop integral function 
C{ti,t2,M^,m^), where one of the three external momenta has zero mass and 
two of the propagators have the same mass M. Here ti = (gi -|- ^2)^, h = and 
q! = 0: 



{27rY (p2 - M2)((p + gi)2 - M^){{p + qi + gs)^ - 

l~x 



dx dy 
Jo Jo 



167r2 Jo Jo M2 - y{A + ti) + xy{ti - ta) + V^h 
Li2 \ — ^-^ I + Li2 



{Annh-t2) { '\y+{t2)) '\y-{t2 
-Li2\ — \—\-Li2 



.y+iii)) \y-ih)J J 

y±{t) = ^{t + A±^A(t,M2,m2)| , (A.14) 

where Li2 is the dilogarithm 

Lt2{x) = - f'^\og{l-xy) . (A.15) 
Jo y 



B Decomposition of the hadronic tensors I^^ 

Here we consider the tensors 

JM- = J dxe^^^+^'^y (0 I TV:^{x)r,_,,{y) \ K^p)) , I = V,A (B.l) 
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and detail its connection with the matrix element (IST 



The general decomposition of A^'^, V^^ in terms of Lorentz- invariant ampli- 
tudes reads [O, for 7^ 0: 



1 \ {2W^ + q^)W'^ ; 



(Ml- - 1^2)^2 



+ l^rfjf^^Sv^ + ^3) (gW^g^ - q'w^)^ (B.2) 



and 



-Ll/^- = iV.e^'^'^q^pp , (B.3) 

where the form factors H^^) and Vi(g^, H^^) are analytic functions of g^ and 

ly^; Fy{q^) denotes the electromagnetic form factor of the kaon {Fy{{]) = 1). 
The amplitudes I^'^ satisfy the Ward identities (cf. Appendix 0) 

q^A^'' = -V2FKP'' 

q^Vf^" = . (B.4) 

The decomposition in Eq. ( |B.2|) may be derived as follows . One first isolates 
the contribution from Bremsstrahlung off the kaon, 

j^fj.f _ j^fj.i' _|_ ^Mi' 

Af,^ = -V2FK ^f^f^, (2iy^ + q^W^ (B.5) 

q,A^ = -V^FKF^iq^)^ . 

The general structure of the remainder A^^ is 

- A^"" = AiWq^ + AagV + Asq^WW + A^g"" + A^q^Wr (B.6) 

The factor g^, which multiplies ^3, is dictated by the Ward identity Eq. ( [B.4| ) 
which implies 

qWAi + q^A2 + Ai = V2Fk 

q'qWAs + q'A, = V2Fk [l - F^ (q^)] . (B.7) 

Eliminating A4 and A^ gives (p.2| ). 

In the process ( |3.1|) the photon is real. As a consequence, only the two form 
factors Ailo,W'^) and Vi{0,W^) contribute. We set 

A{W^) = Ai{0,W^) 

V{W'^) = Vi{0,W'^) (B.8) 
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and obtain, for the matrix element ( p.2| ): 



(B.9) 



q''=0 



with 



F + miu{p,){l + 



2pf+ 



Mpi) 



(B.IO) 



Grouping terms into an IB and a SD piece gives (|3.2|) , (|3.3|) . As a consequence of 
( |B.4| ), T is invariant under the gauge transformation ^ + q^. 

The amphtudes Ai, A2 and Vi are related to the corresponding quantities 
Fa, i? and Fy used by the PDG by 



:b.ii^ 



The last term in ( p.2| ) is omitted in [jTOl . It contributes to processes with a virtual 
photon, l^uil'+l'-. 



Finally, the relation to the notation used in [11, 12 



IS 



2{A±Vy = {a,±VkY [11] 
V2iA,V) = iFA,Fv) il . 



:b.i2) 



C Ward identities 

Here we derive the Ward identity 

A^''{q,p) = I dxe^'^-^^'^y{Q\TV!;^{x)Ai_,M\K^{p)) (c.i: 



g^A^^ = -V2Fkp'' -^W = p-q (C.2) 
without using formal manipulations with T-products and with current commu- 



tators. Our treatment is similar to the considerations in Ref. and parallels 



the ones in Ref. 35 



The reason we wish to discuss this issue is the following. Owing to the short- 
distance singularities of the operator product of two currents, integrals such as 
the one in Eq. ( |U.1| ) are a priori not well defined. As a result of this, the Ward 
identities for A^^ depend on the convention chosen for ({)\TV^{x)A^{y)\K) . In 
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CHPT, the Green functions are obtained through (functional derivatives of) the 
generating functional Z{v,a, s,p). This functional is constructed such that it is 
invariant under transformations generated by the vector currents (see Section Q), 

Z{v' ,a' , s' ,p ) = Z{v,a, s,p) , (C.3) 

where 

s' + ip' = g{s + ip)g^ 

g{x) e ^f/(3) . (C.4) 

The relation ( |C.3| ) embodies all the Ward identities associated with the trans- 
formations ( p.4| ). In particular, seagulls and Schwinger terms are automatically 
taken into account [01, |35|1. For an infinitesimal transformation 



g = 1 + ia + 0(q;^) , 
61 = i[a,I] ■,I = a^,s,p, (C.5) 



one obtains from (|C.3|) 



where 



The Ward identities for current correlation functions are obtained from this gen- 
eral formula by taking additional functional derivatives, evaluated at 

Vf, = = p = , s = diag(m„, ma, rris) • (C.8) 

In the following we consider the isospin symmetry limit rriu = ma. In order to 
derive Eq. ( |C.2| ), we consider the quantity 

T^''^{x,y,z) = {0\TVi'^{x)Al_,,{y)A^,^,,{z)\0) , (C.9) 

which provides an off-shell extension of A^" , and translate the Ward identity for 
Y/j-i^p i]2to the corresponding relation for A'^'^. Differentiating ( p.6|) with respect 
to al~^^{y) and ap'^^{z) at ( p.8| ), one obtains with a = |diag(2, —1, —1) 

^:T'^''^{x,y,z) = -[6\x~y)-6\x-zmTAUM<+^,im ■ (CIO) 
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Next we introduce the Fourier transforms 



x,y,z) 



Ml - 
Ml - 



where 

By use of [] 

one finds from ( |(J.10D 



K 



hm ■ 



RY = 2Flp''pP 
R^^p = ^FkAp^pP 



(C.ll) 

(C.12) 
(C.13) 
(C.14) 



The other Ward identities considered in this article may be derived in an 
analogous manner. 



^Again note that we work with the Condon-Shortley-de Swart phase convention. 
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